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, Abstract: We study properties of the semilinear elliptic equation Ait = i on domains 

' in R", with an eye toward nonnegative singular solutions as limits of positive smooth 

£NJ . solutions. We prove the nonexistence of such solutions in low dimensions when we also re- 
quire them to be stable for the corresponding variational problem. The problem of finding 

0^ ■ singular solutions is related to the general study of singularities of minimal hypersurfaces 

, of Euclidean space. 

■ 1. Introduction 

£h One way to obtain singular minimal hypersurfaces with symmetry is given 

'— in the paper [14] by Simon. Given positive u : C R n — ► R solving the 
. equation 

> : 

9: « Mu:=j:D i (-=22L=)=—=^= 

h \y/l + \Du\*J Uy/l + \Du\* 

and an m-dimensional closed subgroup T of the orthogonal group in R^, 
with G p := {g(p) ■ g G T} an orbit of maximal volume over p £ S^ -1 , the 
! "symmetric graph" G(u) defined by 

G(u) = {(x,u(x)lo) :xett,iueG p }c R n+N 

& : 

will be stationary with respect to n + m-dimensional volume, and will have 
the same regularity as the function u. Positive solutions u > are smooth. 
^ , However, if a sequence Uj of such solutions converges continuously to a weak 

| solution u > 0, then u will be singular exactly at the points where it is zero. 

The corresponding G(u) will be a singular minimal submanifold. A degree 
theoretic program for obtaining such sequences of solutions is outlined in 
Simon's paper. For more on equation (1), see the survey paper by Dierkes [6]. 
Our goal is to apply a similar program to the equation An = -. The basic 
degree argument is presented in Section 4. 

Notice that if we linearize the left hand side of (1) in the form 



v - >. - DiuDjuDiDjU m 

2^ DiDiU - ^ 1 + \ Du \2 = -> 

the resulting equation is An = ^. For this equation the constant m scales 

with the independent variable, i.e. u(Cx) solves An = ^f^, so we may 

l 



restrict our attention to the case m = 1: 

(2) Au = - 

u 

Both equations (2) and (1), with m = 1, have the particular solution 
u(x) = C\x\ with C = l/\/n — 1. By the results of section 3, this solution 
is indeed the limit of a sequence of positive smooth solutions. 

We note that equation (2) is the Euler-Lagrange equation for the varia- 
tional integral 

(3) F{u) = I l\Du\ 2 + logu 

Jn 2 

while equation (1) has variational integral 

/ u m ^Jl + \Du\ 2 , 
Jn 

which is the n + m-dimensional volume of G(u). Notice that for arbitrary 
positive u, the integral (3) is not bounded below for any positive boundary 
data (p on dVt. Indeed the function u e := tp + £(e — (p), where £ is a smooth 
cutoff function on f2, satisfies J-{u e ) — > — oo as e —> 0. 

Equation (2) also arises in relation to chemical catalyst kinetics (See [1] 
and [5]). Here it is a special case of the more general equation 

(4) Au = u~ a , a>0 

For a 1, (4) is stationary for the variational integral 

(5) / -\Du\ 2 + —^—u 1 ^ 



^2' 1-a 

Note that when < a < 1, this integral is bounded below, and by lower 
semicontinuity, as in [11], minimizers among nonnegative functions with 
given boundary data exist. In fact, most results concerning (4), as in [5], 
are limited to the case < a < 1 and are results on minimizers of the 
variational problem. In this way, (2) is an interesting limiting case. 

An important difference in our method is that the singular limit must 
(weakly) satisfy the PDE Au = - on the whole domain Q. Minimizers of 

(5) need not satisfy the Euler-Lagrange equation (4) on the whole domain, 
in fact they need only solve the free boundary problem on the set {u > 0}. 

Despite the nonexistence of minimizers, we may consider solutions u of 

(2) which are also stable for J-(u) in the sense that F(u + tQ > 

at t=o 

for test functions £• In this case, u satisfies the stability inequality 

(6) / C ~ 2 < [ \D(\ 2 

Jn u 2 J n 

for all test functions £. We call such u "stable solutions." 

In Sections 7, 8, and 9, we prove the main results on stable solutions, 
namely, the Holder continuity of stable solutions, the nonexistence of sin- 
gular stable solutions in dimension less than seven, and an estimate on the 
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size of the singular set of a stable solution. The existence and uniqueness 
of stable solutions is presented in Sections 5 and 6. Here we state the main 
results. 

Theorem 1. Let C M. n and suppose there is a subsolution v with Av > 1/v 
on f2 with boundary values ipo on dQ. Then for any ip > ipo, there is a unique 
stable solution u of Au = l/u on f2 with boundary values <p. 

Theorem 2. For a stable solution u with boundary data p < M , and for 
every < a < 1 and £1 CC 0, there is a constant C(n, M,a,(l) such that 
for all x, y £ £1, \u(x) — u(y)\ < C\x — y\ a . 

Theorem 3. Let 2 < n < 6 and let u be a positive stable solution of Au = - 
on the C 1 ' 1 domain Q, with boundary data (p € C 2 ' a (Q,), \<p\2,a < M, and 
tp > e > 0. Then there is a constant 5 = <5(fi, M, e) such that u> 5 onQ,. 

Corollary 4. In dimension less than seven, there are no singular stable 
solutions of (2). 

Theorem 5. Suppose u is a limit of positive stable solutions of Au = - on 
a domain with singular set A = {u = 0}. Then the Hausdorff dimension 
dim n (A) < n - 4 - 2^2. 

2. Basic Facts 

Positive solutions to equation (2) are subharmonic, so the maximum 
principle implies that they achieve their maximum on the boundary of O. 
However, the difference w = u — v of two solutions satisfies the equation 
Aw + — = 0, and so the maximum principle does not guarantee that a 
positive solution is a unique solution of the Dirichlet problem for its bound- 
ary data. Indeed, solutions of the analogous ordinary differential equation 
u" = - need not be unique, and positive radially symmetric solutions in low 
dimension may not be unique. 

Solutions are invariant under homothetic scaling of the graph. That is, 
if u{x) is a solution on £1, then u(Cx)/C is also a solution on the domain 
0,/C for C > 0. In particular, the conical solution scales under homothety 
of the graph to itself, when centered at the origin. 

Nonnegative limits of positive solutions are singular, as in the following 

Lemma 6. If u > is a weak solution of Au = - in a neighborhood of xq 
with u(xq) = 0, then u is not differ entiable at xq. 

Proof. Suppose we have a weak solution u with u > 0, u(xq) = 0, and u 
differentiable at xq. Then Du(xq) = 0, and for any e > 0, in a small enough 
ball B = Bs(xq) we have < u < e\x — xq\. The weak equation in this ball 
is 

(7) f uAC = [ £ 

Jb Jb u 

where £ is C 2 with compact support in the ball. 
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We choose £ radially symmetric such that < £ < 1, £ = 1 on Bj^] arid 
|AC| < £. Then 



fl « Jb s/2 e|x| (n-l)e \2 
where cj n is the volume of the unit n-ball. But also 

- = / nAC < / e\x\-pr = — \ 

Thus, e 2 > (^rr) 2 g , a contradiction for e small enough. 

Using similar methods, we can derive some basic positivity results. 

Lemma 7. Suppose u is a positive and smooth subsolution, i.e. Au > -, 
on a ball B2 p . Then the following hold: 

(1) \ ! u 2 > u n p n 



P 2 Jb 2p \b p 
B 2p V2 n - 1 

Notice that the second property says that there do not exist solutions on 
any ball that are uniformly small, and any solution defined on all of M n must 
be unbounded. 

3. Asymptotically Conical and Radial Solutions 

In this section it will be convenient to consider the rescaled equation 

n — 1 



Au 



u 

which has the conical solution u(x) = \x\. We will also use the radial variable 
r = \x\. 

Caffarelli, Hardt, and Simon proved in [2] the existence of minimal sur- 
faces asymptotic to minimal cones. Following their proof we get a similar 
result showing the existence of a wide variety of singular solutions of (2) on 
the ball which are asymptotic to our conical solution at the origin. 

The wide variety of singular solutions comes from the existence of solu- 
tions with boundary data which are small perturbations of constant bound- 
ary data equal to one on the sphere dB±. As in [2], we are only able to 
specify the perturbed boundary data of the asymptotic solution in the or- 
thogonal complement of a finite dimensional subspace of L 2 {dB\), which 
depends on the rate at which our solutions will be asymptotic to the cone. 
In this case, the finite dimensional subspace is the span of the first J eigen- 
vectors of the operator —A — (n — 1) on the sphere, where the remaining 

eigenvectors have eigenvalues fi large enough that 1 — ^ + y ^ ( TO ~ 2 ) _|_ ^ > m _ 
The projection onto the orthogonal complement is denoted II j. The result 
may be expressed in terms of the scaled Holder norm on annuli, defined by 
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/=0 r<M<2r x=£y \ x ~ V\ 

r < |y| < 2r 

Theorem 8. Given m > 1 and < a < 1, i/iere exisi e and C depending 
on m,n, and a so that for any function g on dB\ with \g\ C 2, a < e, there 
exists a solution u of An = on B\ with Tlj(u — 1) = lijg on dB\ and 
satisfying for < r < 1/2 

r~ m \u - \x\\ 2a . r < C\g\ C 2, a . 

The proof of this theorem is essentially identical to the argument given 
in [2]. ' 

We call smooth positive solutions of An = 2— which are radially sym- 
metric, i.e. u(x) = n(r), r = \x\, "radial solutions." The resulting ordinary 
differential equation satisfied by n is u rr + ^^^r-Ur — = 0, with the partic- 
ular solution u = r. The following two theorems on radial solutions are the 
most useful for further analysis of the PDE. Theorem 9 is due to Brauner- 
Nicolaenko [1], using bifurcation theory in the context of equation (4). We 
have alternative proofs of these and more general facts using basic ODE 
techniques. 

Theorem 9. For any e > 0, solutions of the ODE problem u rr + ^^u r — 
= with n(0) = e and u'(0) = exist uniquely on [0, oo). These 
solutions satisfy n(r) — r = 0(1), and consequently as e — ► 0, the solutions 
u(r) — ► r uniformly on compact subsets. 

Theorem 10. There exist constants C\ and C% depending on n such that 
on the ball Bi(Q), the Dirichlet problem 

An = 2=± onBx 
u = C on dB\ 

has a solution for C > C\, and has a unique solution for C > Ci- For 
n>7, C 1 = C 2 = 1. 

Thus, the radial conic solution n = |x| is indeed a limit of positive smooth 
solutions. An interesting further result is that these conic solutions are 
stable for n > 7 and unstable 2 < n < 6. We state the result for the original 
equation An = 1/n. 

Lemma 11. The conical solutions u = , , are stable for n > 7 and are 

yjn—l J — 

unstable for 2 < n < 6. 

Proof. This follows from the Hardy inequality with best constant 

(n-2) 2 f C 2 



< [ \Dtf 



4 Jq \x\ 2 

for all C G See [10], or for a simple proof see [7]. 
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4. Degree Construction 



We will use the Leray-Schauder degree with several different setups. For 
the basic theory of the degree on Banach spaces, see [4]. In general we will 
use the Banach space B = C 2 ' a (Cl) and open set U = {u £ B : u > g, \u\2, a < 
Ms} where g is a fixed positive bounded function with positive minimum 5. 
A typical operator T : [0, 1] x U — > B will be defined by Tt(u) = v, where v 
is the solution of 



and ipt are boundary data continuous in t with g < ipt < Ms. We use the 
notation deg(/ — Tt,U,0) for the Leray-Schauder degree invariant for fixed 
points of T t . 

In the following results, all solutions are assumed to be positive. Lemma 12 
comes from a basic Schauder estimate. 

Lemma 12. For each < 5 < 1, M$ can be chosen such that any solution 
u of Au = i with 6 < u < \ satisfies \u\2, a < M$. 

Lemma 13. If<po> C2(n) is constant boundary data on the unit ball in ]R n , 
C2(n) as in Theorem 10, and ifU is convex containing ipo and the solution 
of (2) with data ipo, then deg(I — Tq,U, 0) = 1. 

Proof of Lemma 13: 

By Theorem 10, the radial solution with u = ipo on dB\ is unique. We let 
T t {u) = v be the solution to the problem 



so that the unique solution u above is the unique fixed point of To. From 
u we can scale to u = u(y/l — tx) which satisfies Au = and is unique 
relative to its boundary data, u = u(y/l — t) < <po on dB\. We may then 
geometrically scale u to get u uniquely solving Au = and u = ipo 
on dB\. Note that u > u on B±. Thus, Tt has a unique fixed point for 
all t, and in our Leray-Schauder degree setup, there are no fixed points 
of Tt on the boundary of U for any convex ti containing u and ip$. So, 



deg(J - T ,W,0) = deg(/ - T lt U,0). But T x = C. So, deg(J - T lt U,0) = 
deg(J-C,W,0) = deg(/,W,C) = 1. 



Lemma 14. There exists e(O) such that if <pi < e is boundary data on a 
domain Q, in W 1 , then no solution with boundary data 4>\ exists and deg(I — 



Proof. Lemma 14 follows easily from part 2 of Lemma 7 

Lemma 15. If 0, is an arbitrary domain and the function g is chosen to be 
the maximum of a finite collection of subsolutions, i.e. U = {u G C 2 ' a : u > 
gi, . . . , u > tjffc, \u\2, a < Ms} with Agj. > j^, and if ipi > maxg% is boundary 
data on Q, then deg(I — Tx,U,0) = 1. 





on Bi 
on OB 



i 



T u U,0) =0. 
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Proof: Consider again the map T t {u) = v where v is the solution of 

f Av = $ on Q 
\ v = <fio on dQ 

Since To is constant, it has degree one. Suppose for some t that Tt has a 
fixed point u in dlA. Then An = ^, and for some gj, we have u > gj and 
u(xo) = gj{xo) for some xq. Then, 

A (it - gj) < t9j - U < u - gj > on 0. 

So u — <7j has a zero minimum, contradicting the Hopf Maximum Principle. 
Thus, deg(J- T u U,0) = 1. 

We now outline a general method for producing "singular sequences" of 
positive solutions to (2) with minimum tending to zero. In the application 
of the degree, let us choose = B\ C W 1 . Let g = 5j > so that U = {u E 
C 2,Q : n > 5j, \u\2, a < Mg.}, and let Sj \ 0. Using Lemmas 13 and 14, we 
may take (ft to be any homotopy of boundary data between (fo = C > C(n) 
a large constant and ipo < e small. Then since deg(/ — To,b(,0) = 1 and 
deg(7 — T\, U, 0) = 0, there must exist tj £ (0, 1) and a fixed point Uj £ dU 
which solves 

f Am = — on Si 
[ Uj = (p tj on dBi 
with min Uj = 5j . Then the sequence Uj is a "singular sequence" in the sense 

that miniij — > 0. If Uj — > u uniformly with u > and minu = 0, then u is a 
singular solution. Notice that if ip t G C 1 is bounded, then we at least have 
a subsequence j' and a to with — > ip to ■ However, we do not yet have the 
necessary continuity estimates on Uj to get a singular limit u. 

In the case n = 2, since the Uj are in particular subharmonic, we can use 
the "log trick" (See Corollary 24 below) to show that 

Jb p |logp| 

uniformly as p — > 0, which is just short of a modulus of continuity estimate 
and also shows the Uj are uniformly of vanishing mean oscillation as in [19]. 
Even in dimension 2, subharmonicity cannot be sufficient for a continuity 
estimate. For example, the functions u t {x) = \x\ e for e > are nonnegative, 
uniformly bounded on B\ C M 2 , and subharmonic, yet do not satisfy any 
continuity estimate. 

5. General Facts about Maximal Solutions 

A maximal solution u of Au = ^ satisfies the property that v < u for any 
other solution v with the same boundary data as u. We show in Lemma 17 
below that whenever a subsolution exists for fixed boundary data, there is 
also a maximal solution with that boundary data. It turns out that the 
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maximal solution is also the unique stable solution. The existence of max- 
imal solutions can be achieved by the usual method of sub/supersolutions 
(see [5]). We give an alternative degree method. 

We will consistently use the Leray-Schauder degree with the operator 
T t (u) = v, where v is the solution of 

f Av = \ 

\ v\on = (ft 

Lemma 16. For a nonempty finite set of positive subsolutions Uj with 
boundary data (p±, there is a solution u with boundary data ip\ such that 
u > Uj for all j. 

Proof. Consider the open set 

U = p| {u G C 2 ' a : u(x) > Uj(x), \u\ 2 , a < M 5 ) 
j 

in C 2 ' a (ft), where all Uj > 5 and M$ is chosen according to Lemma 12. 
Take data (po > ip\ and let <pt be any smooth decreasing homotopy from 
ipo to Lp\. By Lemma 15, deg(I — T t ,U,0) = 1 for all t < 1, and thus 
there exist solutions ut for all t < 1. Consider any sequence tj /* 1, 
and corresponding solutions u t . . Since all of these functions are uniformly 
bounded below, the Schauder estimates give us a uniform C 3 bound, so 
by Arzela-Ascoli, a subsequence ut k — > u G U in C2.0, and Am = 1/u, 
u\dn = <Pi, u > uj yj. 

Lemma 17. // there is a positive subsolution uq to Au = - on Q with 
boundary data ipo, and if ip > fo, then there is a unique maximal solution 
with boundary data ip. 

Proof Consider the collection C of all solutions u > uq with boundary data 
ip. By Lemma 15, this collection is nonempty. C is partially ordered by the 
relation u a < up on O. By the Hausdorff Maximality Theorem, there exists 
a maximal totally ordered subset 5. For any xo G O let u a ,up G 5 with 
Ua 7^ up and u a < up. By the maximum principle, we have u a (xo) < up{xo). 
Thus, 5 can be indexed by u a (xo). That is, S = {u a } aGA where a = u a (xo). 
By the maximum principle, since Au a > 0, u a < sup(p for all a, so A is 
bounded above. Let = sup A a. We claim that G A, and u aoo is a 
maximal solution. Consider u aj G S with oij f Ooo and u aj < u aj+1 . Since 
these are uniformly bounded below and monotone increasing, the Schauder 
estimates and Arzela-Ascoli give a function u aao with u a . /* u aoc , where 
u aoo( x o) = «oo and u aoo is a solution. For any u a G S, choose j large 
so that ctj > a. Then u aj > u a by total ordering, and u aac > u a . since 
the sequence was monotone. Thus, u aoc is an upper bound for S, so by 
maximality u aac G S. To see that u aoc is a maximal solution, suppose v is 
another solution with v(x\) > u aoc (x\) for some x\ G Q. By Lemma 16, 
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there exists a solution v with v > u aoo and v > v. Also by the hypothesis on 
x\, v 7^ u Qoo . But then SU {v} is totally ordered, contradicting maximality. 



Lemma 18. If po < <fii> and uq,u\ are maximal solutions with boundary 
data ipo and fi respectively, then uq < u\ . 

Proof: By Lemma 15, there exists a solution u to Au = - with data ip\ and 
u > uq, since uo is a subsolution. Since u\ is maximal, u\ > u > uo- 

Lemma 19. If u is a maximal solution on Q and £1 C 0, is a subdomain 
with continuous boundary, then u restricted to 0, is a maximal solution with 
respect to its boundary data on dti. 

Proof. Suppose not. Then there is a maximal solution v on Q with boundary 
data u and v > u on f2. Let U be the open set 



U 



w £ C 2 ' a : w > u on £1, w > v on f2, |u>|2,a < Mgj 



Let </?i be any boundary data on d£l greater than the boundary data (fo of 
u. As in the proof of Lemma 15, consider the operator T t (w) = w where w 
is the solution of 

Aw = l A on Q 
w = (pi on d£l 

Suppose w is a fixed point of Tt in U. By the Hopf Maximum Principle, 
w > a on n. Thus, tu > u on d£l. Then we may apply the maximum 
principle on f2, so w > v on f2. Thus, w cannot be on the boundary of IA. 
So, deg(I — Ti,U,0) = deg(I — Tq,U,0) = 1. Then, for any boundary data 
(ft > ipo, there exists a solution w of Aw = — on with w = ft on dtt, w > u 
on 17, and w > v on f2. Now we let (ft be any smooth decreasing homotopy 
of boundary data approaching tp^. Let wt be the corresponding solutions 
whose existence we just proved. By the Schauder estimates and Arzela- 
Ascoli, there exists a sequence wt j with tj — > such that Wt converges to a 
solution w with boundary data ipo, and with w > u on f2 and w > i; on O. 
Thus, w > v > -u on f2, contradicting the maximality of u. 

6. Stability of Maximal Solutions 
Recall that stable solutions satisfy the stability inequality (6). 
Lemma 20. Maximal solutions of Au = 1/u are stable. 

Proof. Let uq be a maximal solution with data tpo, and let ipt = <po + t for 
t > 0. By Lemma 17, there exist maximal solutions ut with data ft and 
ut > uq. By the Schauder estimates the ut are also bounded in C 4 . For a 
sequence tj \ 0, we then have a subsequence such that ut J — ► uq i n C 2,a , 



with uq > uq, and by maximality uq = uq. Let Sj = maxfi(u tj —Uq), and let 



vj = — 2 g-, — . Then 

Avj = -i. 

UQU tj dj 

By the Schauder estimates, Vj is bounded in C 2,Q , so by Arzela-Ascoli, a 
subsequence Vj — > v in C 2 . The function v is nonnegative, not identically 
0, has nonnegative boundary data, and satisfies the linearized equation Av + 
■\v = 0. By the maximum principle, v > in $7. The weak equation for u 
is then 

Dv-D(= / V 

Jn Uq 

We use the test function ^— for Q, and Cauchy-Schwartz to get the desired 
inequality 

-2 



for all compactly supported £. 

In fact, the maximal solution for given boundary data is the only stable 
solution. 

Lemma 21. The maximal solution for given boundary data (p is the unique 
stable solution with data (p. 

Proof: Let u be the maximal solution and let v be any other positive solution. 
Then u = v + w where w > in $7 and w = on 90. Thus, 

11 — w —w 
Aw = = — > -^r- 

v + w v v(v + w) v z 

and so, integrating by parts with w, 

,2 



— =- > / wAw = I \Dw 

a v 2 J n 



2 



and v cannot be stable. 

Theorem 1 now follows from Lemmas 17, 20, and 21 



7. Holder Continuity of Stable Solutions 

We now prove a main result that stable solutions are locally uniformly 
Holder continuous. 

Theorem 22. For a stable solution u with boundary data ip < M , and for 

every < a < 1 and (l CC there is a constant C(n, M,a,(l) such that 
for all x, y G O, \u(x) — u(y)\ < C\x — y\ a . 
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Proof: Let a be a smooth positive stable solution on 0. The weak form of 
the equation is 

(8) / Du-D( = - 



u 



for C compactly supported. Here and throughout the rest of the proof, all 
integrals are taken over the domain f2. Substituting u( 2 for ( in (8) yields 



c 2 



J \Du\ 2 ( 2 + 2 J u(Du ■ DC, = - J 
and applying the Cauchy-Schwartz inequality, we get 

(9) |c 2 Ql^| 2 + l)< 2 | U 2 | J DC| 2 . 
Substituting in (8) and again using Cauchy-Schwartz gives 

(10) /^C 2 <2/g + 4/|£C| 2 . 

Differentiating the equation with respect to xi, and using subscripts to de- 
note differentiation, we have Aui = — and thus the weak equation 



-C 



where we sum on the repeated index i. Note that this equation is equivalent 
to using Q in (8) and integrating by parts. Substituting ui( 2 for £ gives 



uuuiiC + 2 / muiCCi = I ^C 2 



or, after summing on I, 



' \D 2 u\ 2 c 2 + 2 / u lUli cCi = I ~^C 2 



\Du 



2 



It is now convenient to use the variable v = Wl + \Du^, where 

I 1 2 ^ i 1 2 



V 

j 



and 

\Dv\ 2 < ^2\uij\ 2 = \D 2 u\ 2 . 

Also, Wi = uiuu. Replacing in equation (11) gives 

\Du 



2 



(12) / \Dv\ 2 C 2 + 2 / vviCCi < j L ~^C 2 

If u is stable, it additionally satisfies (6), 
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Now we can combine our inequalities: 

f\Dv\ 2 e + 2fvv i (( i <J l -^fC 2 by (12) 

<2f£ + 4f\D(\ 2 by(10) 

<6/|DC| 2 by (6) 
and we get the main inequality for stable solutions: 

(13) f \Dv\ 2 Q 2 + 2 j vv t CCi < 6 J \DQ\ 2 

We use this and the Sobolev Inequality to iteratively estimate integrals 
f v q (P for q > 0. In fact, the estimate for q = 2 and f3 = 2 is contained 
in (9). For q > 2, replace £ in (13) by v 9 C to get 



J\Dv\ 2 v\ 2 + 2qJ\Dv\ 2 v 2 i( 2 + 2j 
< 6 J \qv q ~ l C,Dv + v q DC 



2 



and so 



(14) j \Dv\ 2 v 2q Q 2 + 2q J \Dv\ 2 v 2q C 2 < 

2 J |Z>u||.DC|v 2<r+1 C + I2q 2 J \Dv\ 2 v 2g - 2 ( 2 + l2 J v 2q \D(\ 2 

where we have used the squared triangle inequality (a + b) 2 < 2a 2 + 2b 2 . We 
use the Cauchy-Schwartz inequality to eliminate the second term on the top 
line of (14) with the first term on the bottom line. 



(15) 



J\Dv\ 2 v 2q e < 1- J V 2q + 2 \D(\ 2 + 12 J v 2q \D(\ 2 + 
+ I2q 2 J \Dv\ 2 v 2q ~ 2 C 2 



Notice at this point that the last term on the right side of the inequality is 
the same as the left hand side with a lower power of v. So, we can apply (15) 
to that term iteratively until 2q — 2 is less than zero, and use the fact that 
v > 1 to get 

(16) j \Dv\ 2 v 2q C 2 < C(n,q) J v 2q+2 \D(\ 2 . 
We rewrite the equation above as 

(17) j \D(v q+1 ()\ 2 < C(n,q) J v 2q+2 \DC,\ 2 , 
replace q by q — 1 , and apply the Sobolev inequality to get 

(18) (J v 2qK ( 2 A K < C(n,q) J v 2q \D(\ 2 for q > 1 
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with At = ^2 or k = 2 if n = 2. We now replace £ by and we fix (" so 
that |-DC| 2 is bounded pointwise by C(Cl) and ( = 1 on O. 



(19) fj v 2 q ^\ 2 - < ^ (J v 2 H 2 P -2\ 29 

where the constant C now depends on n, q, £1, and (3. Now we iterate the 
inequality (19) with q = 1,K, n 2 , . . . and corresponding [3 = 2, /3i,/?2,---, 
where /?,• = 1 + k + k 2 + • • • + 2k>' . Then we have 



(20) n>«» C 2(2K»+K"- i +...+«A 2K,n < c(n)m) o) (y x 5 

By (9), 



i i 

< C{n,m,n,M) 
So, on the subdomain we now have a bound for the Sobolev norm: 
(21) \\u\\ wl>2K m^ <C(n,m,0,,M) 

By the Sobolev Imbedding Theorem, to each a in the statement of the 
theorem there corresponds an m in equation (21) depending on n and a 
such that we have a bound on \u\co,a(Qy This completes the proof. 

Remark: This theorem may be just short of a sharp interior regularity 
estimate, since the known conical example solutions are at worst Lipschitz. 
For equation (4) with < a < 1, a sharp estimate for solutions of the free 
boundary problem minimizing the variational integral was given by Phillips 
in [12]. 

8. Lower Bounds for Stable Solutions in Low Dimensions 

Recall that for n > 7, the radial solutions are unique for their Dirichlet 
boundary data, therefore maximal and stable. In particular, the conical 
solution is a stable singular solution for n > 7. However, for 2 < n < 6, the 
conical solution does not satisfy the stability inequality (6) by Lemma 11 
and is not maximal for its boundary data. Thus it cannot be the limit of 
stable radial solutions and, for 2 < n < 6, the stable radial solutions are 
bounded below by a constant. The next result generalizes this lower bound 
to all stable solutions in a compact subdomain. Theorem 25 gives the same 
result on the entire domain. 

Theorem 23. Let 2 < n < 6 and let u be a positive stable solution of 
Au = ^ on the domain Q, with u < M, and let £1 CC be a compact 
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subdomain. Then there is a constant 5 = 5(n, Cl, M) > such that u> 5 on 

n. 

Proof. The proof follows from the estimate 

(22) I u~ p < C(n,p) for p < 4 + 2y/2. 
Jq 

Notice that the restriction on p allows for p > n as long as n < 6. For 
u > smooth, we use the stability inequality (6) with the test function 
(,u~ q . Then for e > 0, 

J u- 2q - 2 C 2 < J \u~ q D( - qu-t-^Dul 2 

< J u~ 2q \DC,\ 2 + 2\q\u- 2q ~ 1 C\Du\\DC\+ q 2 u- 2q - 2 (: 2 \Du\ 2 

(23) <U + W\ J u ~ 2q \D(\ 2 + {q 2 + 2\q\e) J u~ 2q ~ 2 C 2 \Du\ 2 

Again, all integrals in this proof are taken over the domain $7. Notice that 
in every integral the integrand has compact support in f2. We will also use 
the weak form of the equation (8) with the test function ( 2 u~^, (3 > to 
get 

(24) j u-^C^Du] 2 < J u^ 1 C 2 + 2 J u-P(\Du\\D(\ 
and using Cauchy- Schwartz, for any 5 > 0, 

(25) (/? - 25) J u-P-\ 2 \Du\ 2 < J u-^C 2 + ±J u-? +1 \DC\ 2 

Replacing (3 by 2q + 1 and combining with (23), we get for q > ^- and 
e,<5 > 0, 

Then for 1 — \/2 < q < 1 + \/2 and e and 5 small enough depending on q, 
the coefficient in the last term above is less than one. So, 

(26) J ( 2 u~ 2q - 2 < C{q) j u~ 2q \D(;\ 2 
and now assuming q > 0, we replace C by ( q+1 : 

(27) /(i) 2 ' +2 < OW /(Q 2 '| DC p 
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Now we use Young's inequality in the form 

e a „ a — 1 , 
ab < —a a H — 

with a = and replace g by | — 1. Then for 2 < p < 4 + 2\/2, 



(28) 



So we get equation (22) for any p < 4 + 2\/2. We now recall our continuity 
estimate for stable solutions, that for any a < 1, the Holder norm |u| q < 

C7(fi,a,Af). 

Now let O C CC ri with dist(f2, 317) > so that for any x £ J7, B p (x) C 
fj. For xo G ^, let r = |x — xq\ and suppose u(xo) = e. Then for x € 
B p {x ) C fi, u(a?) < e + C(fi, a, M)r a . Then 



p r n-i dr . 



/V p > / vT v > f (e + Cr Q )" p > nuj n f (e + Cr Q )' 
Jn Jb p (x ) Jb p (x ) Jo 

We may choose a and p large enough that n — 1 — ap < —1 for n < 6. Then 
for e small enough, we have a contradiction of equation (22) with f2 in place 
of £1. This completes the theorem. 

We note that we did not need to use the continuity estimate in the above 

proof. In fact, equation (22) together with the LP estimates gives an estimate 

1 1— — ~ / — 
for u G C ' p (Q) for n < 6 and p < 4+2V2. We use a similar method below 

to get a lower bound on the whole domain. First we present an interesting 

corollary. 

Corollary 24. There are no complete stable solutions of Au = ^ on all of 

R n for 2 < n < 6. 

Proof. Suppose not. From (28), for 2 < n < 6 we have 



f(iy<c(„)J\D ( r 



We choose C equal to one on the ball Br, equal to zero outside B R 2, and 
log N 

log if. 



equal to 2 — j^w 011 ^R 2 \ Br- Then, using the variable r = \x\, we have 



±<C[" I"" ■ dr < c 



>B R u n ~ Jr r n {\ogR) n ~ (logi?)™" 1 

and the result follows letting R — > oo. We thank Neshan Wickramasekera 
for pointing out this trick, which also appeared in reference to the Bernstein 
Theorem in [17]. 

Theorem 25. Let 2 < n < 6 and let u be a positive stable solution of 
Au = - on the C 1 ' 1 domain Q, with boundary data ip 6 C 2,a (Q), \(f\2,a < M, 
and if > e > 0. Then there is a constant 5 = 5(Q, M, e) such that u > 5 on 

n. 
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The theorem follows from the following lemma. 



Lemma 26. Let 2 < n < 6 and let u be a positive stable solution of An = - 
on the domain £1, with Lipschitz boundary data tp, and ip > e > 0. Let 
2 < p < 4 + 2\/2. Then there is a constant C(p, \D<p\) such that 

1 < C|0| 



vP - eP 

Remark: Notice that in Lemma 26 there is no assumption on the smooth- 
ness of the domain. 

Proof: Let e > as in the statement and assume u > is a stable solution. 
Let rj > and first consider the stability inequality with the test function 
C= (<p-u-ri) + : 

(<p- u- T]f 



(29) J - 2 1< J \D( V -u-r,) + f 

So, 

J \D (ip — u — n) + | 2 = J Dip D (ip — u — rj) + — J Du D (ip — u — rj) + 
= J Dp-D(p-u- V ) + + f 



< J\D<p\ 2 + lJ\D{ V -u- V ) + \ 2 + lJl (by(29)) 

So, 

(30) J \D (<p - u - n) + | 2 < C J (1 + |ZV| 2 ) 

Now we use (27) with the same test function, and replace q by p/2 — 1 so 
that for 2 < p < 4 + 2^2, 

<*) /(i^)% cw /(<^)'" 2 ,.,,- u) p. 

Recall from equation (25) that for j3 > 1, 

ju-V( 2 \Du\ 2 <C ju-^e + C ju 2 -?\D(\ 2 . 
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p-2 

So, assuming p > 4, we replace (3 by p — 2 and £ by (y? — u — r/) + 2 and we 
have 

'[(p-u-rj) + Y7l l9 . „ f((v-u-vi)_ x 



We use this with (31) to get 



(tp-u-rj) V f {(ip-u-rj),^ 



+ / P — \D(<p-u)\*< 



{ip-u-rj) \ P 2 /■ / {tp - u - rj) x 



where the constant C now depends also on l-D^I- Now we can apply Young's 
inequality twice. Then 



(ip-U-T]) P (if-U- T]) P 

— H 

up vP- 2 



+ „_ 2 " + \D (if - u) | 2 < C / 1 + \D {if - u) 12 



We then use our estimate (30), and let ij tend to zero. 



(v?-u)_ 



v 



< c\n\ 



Now for ip > 2e, 



so 



1 < (C + 1)|0| 



as required. 



Remark: Notice that for p = 2, we get a stronger result y removing the 
dependence ov C on \Dtp\ and using inequalities (29) and (30). Namely, 

(32) / - 2 <- 2 [ (1 + |^| 2 ) 

where we may assume that tp is merely in W 1,2 (Q). 
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Proof of Theorem 25: The lemma demonstrates the inequality 



(33) J±<C(p,n,e,\D<p\) 

But of course ^ = (Au) p by the equation, and we can apply the L p estimates 
(see [9] 9.14) related to the Calderon-Zygmund Inequality. So, 

< C(p,tt,M,e) 

and by the extended Sobolev Embedding Theorem, 

M c i,i-n/j. < C(p,£l,M,e) 
which in particular implies a uniform Lipschitz bound on u. Then if u 
achieves the value 5 at a point xq £ f2, 

i r i 

> 



u p ~ J (5 + Cr)P 

where r = x — xq, a contradiction of (33) for 5 < 5(p, O, M, e) and p > n. 

Remark: With this lower bound we in fact have complete regularity of 
stable solutions for n < 6. From u £ C 1,Q and thus (by the lower bound) 
1 g (7 1 ' ct ) we can apply Holder estimates to get continuous derivatives of all 
orders on the interior of the domain. 

9. Hausdorff Dimension of Singular Sets of Stable Solutions 

We use Hausdorff dimension as described in [15]. 

Theorem 27. Suppose u is a limit of positive stable solutions of Au = - on 
a domain £1 with singular set A = {u = 0}. Then the Hausdorff dimension 
dim n (A) < n - 4 - 2y/2. 

Proof. We will show for any ball B p of radius p whose closure is contained in 
Q, and any (3 > n-4-2 v / 2, that the Hausdorff Measure HP {A n B p / 2 ) < oo. 
First, for any 5 with < 6 < p/4, we cover A n B p / 2 by cubes Qj of side 
length 25 with disjoint interiors, j = 1, . . . ,N. Let p < 4 + 2\/2. By (22), 
we have 

/ — < K < oo 



UQ 3 uP JB p/2 u p 

with K independent of 5. By Theorem 22, for any < a < 1, we have 
u(x) < C (dist(x, A)) a . Thus, assuming all the Qj intersect A, 

C >[ * > f I -> 



Qj ^ ~ Jqj (dist(x, A)) ap ~ J 0<Xi<2S (xf + • • • x z) ap/2 

> J_f 5 rujnr"- 1 ^ > nu n 
- 2 n J r * P - 2 n (n- ap) 
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Choose a so that [3 = n — ap. Then we have 

Wf (A n B p/2 ) <C'J2^ < C"K < oo 

Qi 

independent of 5, and the result follows. 

Remark: Recall that the analogous equation (4) with < a < 1 has actual 
minimizers for the variational problem. In [13], Phillips proved a Hausdorff 
estimate on the free boundary for minimizers. Note that in the free boundary 
problem, the solution is allowed to vanish completely and not satisfy the 
PDE on an interior set of positive measure. The technique above gives an 
estimate on the size of the singular set of a singular solution which is a limit 
of positive solutions satisfying the PDE on the whole domain. 

Remark: The major results of this paper extend to the equation (4) with 
< a < 1, except for two points. Solutions of (4) which achieve the value 
zero need not be singular, and the results of sections 8 and 9 are more 
complicated, with dimensions depending on a. 
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